Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory.

Polariton chemistry has attracted great attention as a potential route to modify chemical structure, properties, and reactivity through strong interactions among molecular electronic, vibrational, or rovibrational degrees of freedom. A rigorous theoretical treatment of molecular polaritons requires the treatment of matter and photon degrees of freedom on equal quantum mechanical footing. In the limit of molecular electronic strong or ultrastrong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of ab initio quantum chemistry, yielding an approach we refer to as ab initio cavity quantum electrodynamics, where the photon degrees of freedom are treated at the level of cavity quantum electrodynamics. Here, we present an approach called Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction theory to provide ground- and excited-state polaritonic surfaces with a balanced description of strong correlation effects among electronic and photonic degrees of freedom. This method provides a platform for ab initio cavity quantum electrodynamics when both strong electron correlation and strong light-matter coupling are important and is an important step toward computational approaches that yield multiple polaritonic potential energy surfaces and couplings that can be leveraged for ab initio molecular dynamics simulations of polariton chemistry.

G0W0 Ionization Potentials of First-Row Transition Metal Aqua Ions.

We report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V-Cu) aqua ions in their most common oxidation states, i.e., V3+, Cr3+, Cr2+, Mn2+, Fe3+, Fe2+, Co2+, Ni2+, and Cu2+. The d-orbital occupancy of these systems spans a broad range from d2 to d9. All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G0W0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G0W0@PBE, G0W0@PBE0, and G0W0@r2SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G0W0@PBE0 and G0W0@r2SCAN approaches, with the r2SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

PTM-Psi: A python package to facilitate the computational investigation of post-translational modification on protein structures and their impacts on dynamics and functions.

Post-translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python-based workflow, called PTM-Psi, integrates several established open-source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non-standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S-nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM-Psi for interpreting sequence-structure-function relationships derived from thiol redox proteomics data. We demonstrate that the S-nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

NWChem: Recent and Ongoing Developments.

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

Basis Set Selection for Molecular Core-Level GW Calculations.

The GW approximation has been recently gaining popularity among the methods for simulating molecular core-level X-ray photoemission spectra. Traditionally, Gaussian-type orbital GW core-level binding energies have been computed using either the cc-pVnZ or def2-nZVP basis set families, extrapolating the obtained results to the complete basis set limit, followed by an element-specific relativistic correction. Despite achieving rather good accuracy, it has been previously stated that these binding energies are chronically underestimated. In the present work, we show that those previous studies obtained results that were not well-converged with respect to the basis set size and that, once basis set convergence is achieved, there seems to be no such underestimation. Standard techniques known to offer a good cost-accuracy ratio in other theories demonstrate that the cc-pVnZ and def2-nZVP families exhibit contraction errors and might lead to unreliable complete basis set extrapolations for absolute binding energies, often deviating about 200-500 meV from the putative basis set limit found in this work. On the other hand, uncontracted versions of these basis sets offer vastly improved convergence. Even faster convergence can be obtained using core-rich property-optimized basis set families like pcSseg-n, pcJ-n, and ccX-nZ. Finally, we also show that the improvement observed for core properties using these specialized basis sets does not degrade their description of valence excitations: vertical ionization potentials and electron affinities computed with these basis sets converge as fast as the ones obtained with the aug-cc-pVnZ family, thus offering a balanced description of both core and valence regions.

Scalable Molecular GW Calculations: Valence and Core Spectra.

We present a scalable implementation of the GW approximation using Gaussian atomic orbitals to study the valence and core ionization spectroscopies of molecules. The implementation of the standard spectral decomposition approach to the screened-Coulomb interaction, as well as a contour-deformation method, is described. We have implemented both of these approaches using the robust variational fitting approximation to the four-center electron repulsion integrals. We have utilized the MINRES solver with the contour-deformation approach to reduce the computational scaling by 1 order of magnitude. A complex heuristic in the quasiparticle equation solver further allows a speed-up of the computation of core and semicore ionization energies. Benchmark tests using the GW100 and CORE65 data sets and the carbon 1s binding energy of the well-studied ethyl trifluoroacetate, or ESCA molecule, were performed to validate the accuracy of our implementation. We also demonstrate and discuss the parallel performance and computational scaling of our implementation using a range of water clusters of increasing size.

Spin-Crossover from a Well-Behaved, Low-Cost meta-GGA Density Functional.

The recent major modification, r2SCAN, of the SCAN (strongly constrained and appropriately normed) meta-GGA exchange-correlation functional is shown to give substantially better spin-crossover electronic energies (high spin minus low spin) on a benchmark data set than the original SCAN as well as on some Fe complexes. The deorbitalized counterpart r2SCAN-L is almost as good as SCAN and much faster in periodically bounded systems. A combination strategy for the balanced treatment of molecular and periodic spin-crossover therefore is recommended.

Multicomponent density functional theory with density fitting.

Multicomponent Density Functional Theory (MDFT) is a promising methodology to incorporate nuclear quantum effects, such as zero-point energy or tunneling, or to simulate other types of particles such as muons or positrons using particle densities as basic quantities. As for standard electronic DFT, a still ongoing challenge is to achieve the most efficient implementations. We introduce a multicomponent DFT implementation within the framework of auxiliary DFT, focusing on molecular systems comprising electrons and quantum protons. We introduce a dual variational procedure to determine auxiliary electron and proton densities which leads to a succession of approximate energy expressions. Electronic and protonic fitted densities are employed in (i) electron-electron, proton-proton, and electron-proton classical Coulomb interactions and (ii) electron exchange-correlation, proton-proton exchange, and electron-proton correlation (EPC) potentials. If needed, exact exchange among electrons or among protons is computed by the variational fitting of the corresponding Fock potential. The implementation is carried out in deMon2k. We test various electron proton correlation functionals on proton affinities. We find that auxiliary densities can be safely used in electron-electron, proton-proton, and electron-proton classical Coulomb interactions, as well as in EPC, albeit with some precautions related to the choice of the electronic auxiliary basis set that must be flexible enough. Computational tests reported indicate that introduction of density fitting in MDFT is clearly advantageous in terms of computational effort with good scaling properties with respect to the number of electrons and protons treated at the DFT level.

Range-Separated Hybrid Functionals with Variational Fitted Exact Exchange.

This work presents a variationally fitted long-range exact exchange algorithm that can be used for the computation of range-separated hybrid density functionals in the linear combination of Gaussian type orbital (LCGTO) approximation. The obtained LCGTO energy and gradient expressions are free of four-center integrals and employ modified three-center integral recurrence relations to obtain optimal computational performance. The accuracy and performance of selected range-separated hybrid functionals with variational fitted long-range exact exchange are analyzed and discussed. A parallel implementation is presented and benchmarked.

Analytic second derivatives from auxiliary density perturbation theory.

The working equations for the calculation of analytic second energy derivatives in the framework of auxiliary density functional theory (ADFT) are presented. The needed perturbations are calculated with auxiliary density perturbation theory (ADPT) which is extended to perturbation dependent basis and auxiliary functions sets. The obtained ADPT equation systems are solved with the Eirola-Nevanlinna algorithm. The newly developed analytic second ADFT energy derivative approach was implemented in deMon2k and validated with respect to the corresponding finite difference approach by calculating the harmonic frequencies of small molecules. Good agreement between these two methodologies is found. To analyze the scaling of the new analytic second ADFT energy derivatives with respect to the number of processors in parallel runs, the harmonic frequencies of the carbon fullerene C240 are calculated with varying numbers of processors. Fair scaling up to 720 processors was found. As showcase applications, symmetry unrestricted optimization and frequency analyses of icosahedral carbon fullerenes with up to 960 atoms are presented.

Robust and Efficient Auxiliary Density Perturbation Theory Calculations.

A new iterative solver for the recently developed time-dependent auxiliary density perturbation theory is presented. It is based on the Eirola-Nevanlinna algorithm for large nonsymmetric linear equation systems. The new methodology is validated by static and dynamic polarizability calculations of small molecules. Comparison between the analytic and iterative solutions of the response equation system shows excellent agreement for the calculated static and dynamic polarizabilities. The new iterative solver reduces the formal scaling from [Symbol: see text](N(4)) to [Symbol: see text](N(3)). Furthermore, the observed computational scaling for linear alkane chains is N(1.6). This subquadratic behavior is possible in systems with a few hundred atoms because of the very small prefactors of the [Symbol: see text](N(3)) and [Symbol: see text](N(2)) steps remaining in the iterative solver. To demonstrate the potential of this new methodology, static polarizabilities of giant fullerenes up to C960, with more than 14,000 basis functions, are calculated.

QM/MM calculations with deMon2k.

The density functional code deMon2k employs a fitted density throughout (Auxiliary Density Functional Theory), which offers a great speed advantage without sacrificing necessary accuracy. Powerful Quantum Mechanical/Molecular Mechanical (QM/MM) approaches are reviewed. Following an overview of the basic features of deMon2k that make it efficient while retaining accuracy, three QM/MM implementations are compared and contrasted. In the first, deMon2k is interfaced with the CHARMM MM code (CHARMM-deMon2k); in the second MM is coded directly within the deMon2k software; and in the third the Chemistry in Ruby (Cuby) wrapper is used to drive the calculations. Cuby is also used in the context of constrained-DFT/MM calculations. Each of these implementations is described briefly; pros and cons are discussed and a few recent applications are described briefly. Applications include solvated ions and biomolecules, polyglutamine peptides important in polyQ neurodegenerative diseases, copper monooxygenases and ultra-rapid electron transfer in cryptochromes.

Robust and efficient variational fitting of Fock exchange.

We propose a new variational fitting approach for Fock exchange that requires only the calculation of analytical three-center electron repulsion integrals. It relies on localized molecular orbitals and Hermite Gaussian auxiliary functions. The working equations along with a detailed description of the implementation are presented. The computational performance of the new algorithm is analyzed by benchmark calculations on systems with different dimensionality. Comparison with standard four-center and three-center electron repulsion integral Hartree-Fock calculations shows an excellent accuracy-performance relation.